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Abstract 

We derive QCD sum rules for heavy baryons at leading order in 1/mq and at 
next-to-leading order in as- The calculation involves the evaluation of four different 
perturbative three-loop diagrams which determine the Qg-corrections to the Wilson 
coefficients of the leading term in the Operator Product Expansion (OPE). From 
the sum rules we obtain estimates for the masses and the residues of the heavy 
baryons Aq and £q. The perturbative 0(as) corrections to the leading order 
spectral function amount to about 100%, and they shift the calculated values for 
the baryon masses slightly upward. The residues are shifted upward by about 
20 — 50%. For the bound state energy A given by the difference of the heavy baryon 
mass and the pole mass of the heavy quark ttiq we obtain m\ Q — txiq = 780 MeV 
and mx Q - m Q = 950 MeV. For the residues we find \F A \ = 0.028 GeV 3 and 
\F S \ ^ 0.039 GeV 3 . 
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1 Introduction 



There has been a great deal of interest in the physics of heavy hadrons containing one 
heavy quark. The Heavy Quark Effective Theory (HQET) allows one to study the prop- 
erties of the heavy hadrons in a systematic 1/mg-expansion. The leading term of the 
expansion gives rise to the spin-flavour symmetry of Heavy Quark Symmetry (HQS). The 
corrections to the leading HQS results are determined by the small expansion parameter 
Aqcd/^q? where Aq CD 300 MeV is the scale of low energy physics (for a review of 
HQET see for a review of HQS and the sum rule approach for heavy mesons see [0). 

Among the well-known predictions of HQS are e.g. the relations between different 
heavy hadron transition form factors. Take for example, the A& — * A c electro-weak 
transitions. The six form-factors describing this transition are reduced to one univer- 
sal Isgur-Wise function in the HQS limit [§, |], |5|. Even then one still remains with many 
non-perturbative parameters characterizing the process and the heavy baryons partici- 
pating in it. These concern the functional behaviour of the Isgur-Wise function itself, 
the masses and residues of the heavy baryons, and, at next-to-leading order in the heavy 
mass expansion, the average kinetic and chromomagnetic energy of the heavy quark in 
the heavy baryon. All these non-perturbative parameters can be determined by using 
non-perturbative methods as e.g. lattice calculations, QCD sum rule methods || or, in a 
less fundamental approach, by using potential models. 

In the present paper we study the correlator of two heavy baryon currents in the HQS 
limit when uiq — > oo. Using the QCD sum rule method we calculate the masses and 
residues of the heavy baryons associated with the heavy baryon currents. In its original 
form the QCD sum rule method was suggested by Shifman et al. || as a tool to investigate 
the properties of light meson systems. Later on the method was extended to the case of 
light baryons in [0, ||, [10| . The QCD sum rule approach has proven itself to be a very 
powerful non-perturbative QCD-based tool which takes into account the properties of the 
QCD vacuum. It allows one to obtain reliable estimates for hadron masses, their residues 
and their elastic as well as their transition form factors. 

In the heavy-light sector the first leading order analysis (leading both in I/tuq as well 
as in ag) of heavy meson properties within the QCD sum rule approach was performed 
in |[T| . Later on the heavy meson sum rule calculation was extended to include next-to- 
leading order radiative corrections. The next-to-leading order corrections proved to be 
rather important \T2L O, O]. QCD sum rules for baryons with large but finite masses 



rriQ were first studied in fLol , jig] . Later on the methods of HQET were incorporated in 
the sum rule analysis. The leading order QCD sum rules for heavy baryons were first 
considered in |ll|, [L7], |18| , again to leading order both in 1/mg as well as in as- Finite 
mass corrections to these sum rules were discussed in [jl9| . 

In order to improve on the accuracy of the existing QCD sum rule analysis of heavy 
baryons one needs to avail of the next-to-leading order radiative corrections to the sum 
rules. This forms the subject of the present paper. We calculate the QCD radiative 
corrections to the leading perturbative term in the Operator Product Expansion (OPE) 
and, from these, we derive next-to-leading order QCD sum rules for heavy baryons in the 
HQS limit. We then go on to analyze the sum rules and compute the values of the masses 
and the residues of the heavy baryons at next-to-leading order accuracy. 
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The paper is organized as follows. In Sec. 2 we introduce heavy baryon currents 
as interpolating fields for the heavy ground state baryons. In Sec. 3 we construct the 
correlator of two heavy baryon currents by means of the OPE and define the spectral 
density. In Sec. 4 we present our results on the radiative corrections to the perturbative 
part of the spectral density and construct renormalization group invariant QCD sum rules 
by recapitulating some known results on the one- and two-loop anomalous dimensions of 
the currents. Sec. 5 contains the results of our numerical analysis. Sec. 6, finally, contains 
our summary and our conclusions. In Appendix A we provide a detailed collection of 
results on the calculation of the two- and three-loop contributions to the correlator of two 
heavy baryon currents. These results are quite general in that they are given for general 
space-time dimensions and for a general baryonic current structure. 

2 Baryonic currents 

The currents of the heavy baryon Aq and the heavy quark spin baryon doublet {£q, £q} 
are associated with the spin-parity quantum numbers j p = + and j p = 1 + for the light 
diquark system with antisymmetric and symmetric flavour structure, respectively. Adding 
the heavy quark to the light quark system, one obtains j p = | + for the Aq baryon and 
the pair of degenerate states j p = \ + and j p = | + for the baryons Eq and Eg. The 
general structure of the heavy baryon currents has the form (see e.g. |17| and refs. therein) 

J = [q lT CTTq^]T'Q k e ljk . (1) 

Here the index T means transposition, C is the charge conjugation matrix with the 
properties OyJC -1 = — 7 M and C'jJC^ 1 = 75, i, j, k are colour indices and r is a matrix in 
flavour space. The effective static field of the heavy quark is denoted by Q. For each of the 
ground state baryons there are two independent interpolating currents Ji and J2 which 
both have the appropriate quantum numbers to interpolate to the respective ground state 
baryons. They are given by Jy], |l7j 

Jai = [q iT CT^q j }Q k e ijkl J A2 = [q iT CT^ q j }Q k e ijk , 
Jm = [q iT Cr^q J ] ■ 775<2%'fc, </s2 = [q iT Cr^q j ] ■ 775<2%'fc, 
4i = [q lT Cr^}Q k e l3k + ^[q iT Cr^] ■ jQ k e ijk , (2) 

Js*2 = [q iT CT l0 ^q j }Q k e ijk + lp[q iT C^q j } ■ ^Q k e ijk , 

where Js*i and Js*2 satisfy the spin-3/2 condition 7Js*i = (i — 1, 2). The flavour matrix 
t is antisymmetric for Aq and symmetric for the heavy quark spin doublet {Sq, Sq}. The 
currents written down in Eq. (0) are rest frame currents. The corresponding expressions in 
a general frame moving with velocity four-vector v M can be obtained by the substitutions 
7o —>■ 1) and 7 — > 7^ = 7^ — fv^. In the following analysis we shall be using both of these 
equivalent descriptions alternatively, i.e. we shall either use the static description with 

= (1, 0, 0, 0) or a moving frame description with = (l,v) and v 7^ 0. 

For a general analysis it proves to be convenient to represent the general light-side 
Dirac structure of the currents in Eq. (Q) by an antisymmetrized product of n Dirac 
matrices T = 7^ • • -7^. When calculating the one- and two-loop vertex corrections one 
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Table 1: Specific values of the parameter pair 



for particular cases of the light-side 



Dirac structure 17 7^° refers to the naive 75-scheme with an anticommuting 75 0] and 



75 to the 75-scheme due to Breitenlohner, Maison, 't Hooft and Veltman El . 



encounters 7-contractions of the form 'y a T'y a and 7(^70. The 7 a -contraction leads to an 
n-dependence according to 

7 Q r 7 Q = hY = (-l) n (D - 2n)r. (3) 

The 70-contraction depends in addition on an additional parameter s which takes the 
value (s = +1) and (s = —1) for an even or odd number of 7o's in T, respectively. The 
7o- contract ion reads 

7or 7 o = (-l) n sT. (4) 

In order to facilitate the use of Eqs. @ and (|) we have compiled a table of the (n, s)- 
values relevant for the heavy baryon currents treated in this paper (see Table 1). 



3 Correlator of two baryonic currents 

In this section we describe the steps needed for the evaluation of baryonic QCD sum rules. 
One starts with the correlator of two baryonic currents, 

U(u — k-v) —i J d 4 xe ikx (0\T{J(x), J(0)}|0), (5) 

where fc M and are the residual and full momentum of the heavy quark and v^ is the four- 
velocity using the momentum expansion p M = mQV^ + k^. As was mentioned before, there 
are two possible choices of interpolating currents for each of the heavy baryon states, given 
by Ti and T 2 = I\^. Thus one may consider correlators of the same currents (diagonal 
correlators) or of different currents (non-diagonal correlators). In the general case, one 
may even consider correlators built from a linear combination J = Ji + 6J 2 of these 
currents with an arbitrary coefficient b. We mention that the choice 6=1 corresponds 
to a constituent quark model current which has maximal overlap with the ground state 
baryons in the constituent quark model picture. In this paper we limit our attention to 
diagonal correlators only. 

The correlator in Eq. (|5|) depends only on the energy variable uo = k-v because of the 
static nature of the heavy propagator. It can be factorized into a spinor dependent piece 
and a scalar correlator function P(u>) according to 

n(w) = r'i±^F' \ Tr(IT)2 Tr(rrt)P(^). (6) 
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Following the standard QCD sum rule method ||, the correlator is calculated in the region 
— uj ^1 — 2 GeV, including perturbative and non-perturbative contributions, where the 
non-perturbative contributions can in general be quite large. The non-perturbative effects 
are taken into account by employing an Operator Product Expansion (OPE) for the time- 
ordered product of currents in Eq. (|5|). One then has 

T{J(x),J(0)} = ^C d (x 2 )O d , (7) 
d 

where the operators O d are local operators with a given dimension d, Oq = 1, O3 = (qq), 
O4 = (GG), ... , and the expansion coefficients C d (x 2 ) are the corresponding coefficient 
functions or Wilson coefficients of the OPE. 

A straightforward dimensional analysis shows that the OPE of the diagonal correlator 
contains only even-dimensional terms. We take into account the perturbative term for 
d = 0, the gluon condensate term for d = 4 and a condensate term with four quark fields 
for d = 6. The four-quark operator will be factorized into a product of two two-quark 
operators, (q{0)q(x)) 2 |J. Accordingly the Fourier transform of the scalar correlator 
function P(uj) reads 

P(t) = W = «»(W(i + 32 ^f\ V " W! mm ^ <8) 

where c = 1 for Aq, c = —1/3 for {Eq,Eq} and N c is the number of colours. For the 
non-local quark condensate (q(0)q(t)) one may use the OPE about (qq) := (g(O)g(O)), 
namely 



(q(0) q (t)) = (q q )(l + -m 2 t 2 + na s (GG)— + ... ), (9) 

where the parameter m is defined in Eq. (|i~T|). Alternatively one may use the Gaussian 
ansatz [2^ 



(q(0)q(jt)) = (qq)ex V (^mlt 2 ). (10) 

When expanding the Gaussian ansatz one sees that the two forms agree up to the term 
linear in t 2 . Thus the two representations of the non-local quark condensate are quite 
similar to one another for small values of t. In our sum rule analysis we shall make use of 
the Gaussian ansatz because it provides for better stability of the sum rules. 
For the condensates we use the numerical values 



= -(0.23 GeV) 3 , 
a s (GG) = 0.04 GeV 4 , (11) 
9s(q^G^q) = m 2 (qq) with ml = 0.8 GeV 2 . 

With these condensate values one sees that the OPE in Eq. (|8]) with Euclidian time r = it 
converges nicely for 1/r > 0.3 GeV. In this region one may thus safely truncate the OPE 
series after the second term. At 1/r = 0.3 GeV the contribution of the first term is two 
times larger than the last quark condensate term. Its contribution grows rapidly with 
1/r. When 1/r is further increased we see that the correlator becomes dominated by the 
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perturbative contribution. For example, at 1/r = 0.6 GeV the perturbative term is two 
orders of magnitude larger than the contribution of the condensate terms. Note, however, 
that at 1/r = 0.4 GeV the contribution of the ground state to the correlator is ten times 
smaller than the contribution of the excited states and the continuum. This would imply 
that if the theoretical and phenomenological continuum contributions differ by about 10% 
(and are not equal to each other as assumed here), this difference would induce a 100% 
change in the contribution of the ground state. Thus the sum rules can only be trusted at 
values 1/r < 0.4 GeV (see also the discussions of the numerical results of the sum rules). 
In the next section we will show that the perturbative corrections become even more 
important at small Euclidian distances in comparison to the non-perturbative condensate 
contributions. 

As a next step one determines the spectral density using the coordinate space repre- 
sentation P(t) of the current correlator. The simplest way to proceed is as follows. The 
scalar correlator function Pqpe(^) satisfies a dispersion relation 

p 0PE M = />M=f^£_ + PV), d2) 

Jo uj 1 — oo — «0 

where p(u) = Im(P(u))/ir is the spectral density and P'{uj) is a polynomial in u, which 
takes into account possible subtractions in the dispersion representation. The Fourier 
transform of the polynomial P'(uj) consists of the 5-function S(t) and derivatives S^ n \t) 
of the 5-function. A comparison with Eq. (||) shows that one does not in fact need any 
subtractions. We therefore set P'{uj) = 0. Taking the Fourier transform of Eq. (|T2|) 
according to 

P(t) = [ p.e-**P(u>) (13) 

we obtain 

P(t) = i / p(u)e- iuJt du. (14) 



o 

Then we analytically continue P(t) from t > to imaginary times by introducing the 
Euclidian time r = it. After this transformation, Eq. (|TJ) becomes the well known 
Laplace transformation. One may thus use an inverse Laplace transformation in order to 
obtain an Euclidean time representation of the spectral density, 

1 rc+ioo 

p{u) = — / P(-ir)e^dr, (15) 

where c is to be chosen as a real constant to the right of all singularities of P{t). It 
is then easy to check that the form P(t) = 6{t)/t n gives the spectral density piu) = 
z n+1 0(o;)a; n - 1 /(ri - 1), whereas P(t) = 6{t)t n results in p{uo) = -{-i) n - 1 5^ {u) for n > 0. 
Following the argumentation in || we do not include forms of the second kind into the 
spectral density p. So the leading order perturbative contribution and the next-to-leading 
order contribution of the gluon condensate to the spectral density are given by 

p(uo) = p (uj) + p 4 (u), where (16) 

u 5 a s (GG) 
PoH = ^ and Pi ( u ) = c-^- U . 
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4 Radiative correction to the perturbative term 



Next we consider radiative corrections to the leading order spectral density in Eq. (|16D. 
There are altogether four different three-loop graphs that contribute to the correlator 
of two baryonic currents, which are shown in Fig. 1. Contrary to the experience in the 
two- loop case, the most convenient way to calculate the three-loop contributions is to 
evaluate them in momentum space. The fact that all graphs in Fig. 1 have two-point 
two-loop subgraphs greatly simplifies the calculational task. One can first evaluate the 
respective subgraphs such that one remains with a one-loop integration. The subgraph 
two-loop integration can be performed by using the algebraic methods described in [p^j . It 
is important to note that the results of the two-loop integration can be expressed in terms 
of a polynomial function of the external momentum that flows into the subgraph. Hence, 
the remaining integration is a one-loop type integration, where the power of one of the 
propagators has become a non-integer number due to the use of dimensional integration. 
The upshot of this is that all steps of the three-loop integration can be reduced to purely 
algebraic manipulations. 

We present the results of calculating the two-loop and three-loop contributions to the 
correlator in the form 

4 

P(cu) = \ C B + \ 1 J2C l B l , (17) 

i=l 

where we have used the abbreviations A = (-2cu/ / u) (2D_3) , Ai = g 2 s / '(4n) c '(—2a; / '/i) (3jD ~ 7) , 
and where D = 4 — 2e is the space-time dimension. Concerning the colour structure we 
have defined the colour factors Cj (i = 0, . . . , 4) according to the labelling of the graphs 
in Fig. 1. Their values are given by C = N c \, C\ = C 2 = —NJCb and C 3 = C 4 = NJCp, 
where C F = (iV c 2 — 1)/2N C and C B = [N c + 1)/2N C . Values for the scalar coefficients B i 
defined in Eq. ( fP7D are listed in Appendix A. 

Putting everything together, the two-loop and three-loop scalar correlation factor P(u) 
defined in Eq. (|17D is given by 

PM = 



32a; 5 \ f -2w\ 4e 1 (l 107\ 



(4tt 



TT -+ — 



/I J 40 \e 15 J 



a s (n 2 - 4n + 6 40((2) + 61n 2 - 234n + 396 (n - 2)s 

+ \ 45? + 225e + 90 

5(195n 2 - 780n + 1946)C(2) - 2200C(3) + 4907n 2 - 18408n + 34352\ 
+ 2250 ) 

The scalar correlation function P(u) is renormalized by the square of the renormalization 



factor Zj of the baryonic current derived in |jl7 |. Accordingly one has 



P{u) = Z 2 jP rcn {u) with Zj = 1 + ^^(n 2 - 4n + 6). (19) 

The multiplication of P(u>) in Eq. (0) with Z 2 results in the cancellation of the second 
power in 1/e. The remaining 1 / e-singularity is purely real and hence does not contribute 
to the spectral density. Since the renormalized spectral density p ren (u) = Im(P ren (uj))/TT 
has to be finite, this provides a check on our calculation. The spectral density can be 
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read off from Eq. ( |T8"D and is given by 



p ven (uj,p) = p {uj) 



1 + if 



where (20) 



PoM = and r(u/ii) = r x In ( — J + r 2 with 

r x := -(n 2 -4n + 6) and r 2 := — (60C(2) + 38n 2 - 137n + 273). 
3 45 

The coefficient r\ of the logarithmic term in Eq. (p0| ) coincides with twice the one-loop 
anomalous dimension given in Eq. fl36|), as expected. The reason is that the evolution 
of p(u,p) is controlled by the renormalization group equation and that the anomalous 
dimension of (J J) and p(uj,p) coincide. 

The as-correction can be seen to depend on the properties of the light-side Dirac 
matrix T in the heavy baryon current, as specified in Table 1. As an explicit result we 
list representations of the r (oo/fi) -functions of the four baryon currents in the naively 
anticommuting 75-scheme (AC). They read 



r Al (u/p) = 16 In 



p\ 8(20C(2) + 91) 



2u J 15 



S66.04 



«48.48 

8, / p \ 8(60C(2) + 151) 



r^M = gln(£) + 



45 



^44.40 



The results for the two different baryon currents A x and A 2 in the 't Hooft-Veltman 75- 
scheme (HV) differ from those presented above. It is well known that currents in different 
75-schemes are connected by a finite renormalization factor Z such that 

J ac = ZJhv- (22) 
These finite factors also appear in the calculation of two-loop anomalous dimensions of 



baryonic currents |24| . From the results of |24| one has 



Z A1 = 1 - and Z A2 = 1 - (23) 

Using these finite renormalization factors one may convert the results in the naively 
anticommuting 75-scheme given in Eq. (^T|) to the corresponding results in the 't Hooft- 
Veltman scheme. Least the reader worry that we do not list the corresponding S-type 
conversion factors we remind him that the 75 in the S-type currents act on the heavy 
quark line and thus there are no 75-ambiguities. Nevertheless, the 't Hooft-Velman 75- 
scheme needs some counter terms to satisfy some kind of Ward identities. To avoid this 
complification, we will henceforth concentrate on the naively anticommuting 75-scheme, 
where such counterterms are not necessary at all. We only mention that the finite renor- 
malization in Eq. will bring the results of the two IVschemes in line. 
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In order to allow for a quick appraisal of the importance of the perturbative corrections 
we have exhibited the numerical values of the second terms in Eq. ([H]). For as we use 
the running coupling constant, which we normalized to the value of as(mz) = 0.118 at 
the mass of the Z-boson for Nf = 5 active flavours. By doing so one has as{^) = 0.333 
at n = 1 GeV for Nf = 3 active flavours. Using this value for a${^ = 1 GeV), the 
above results show that the perturbative a^-corrections to the spectral density amount to 
about 100%. This highlights the importance of perturbative QCD radiative corrections in 
QCD sum rule applications. The same observation was made in the heavy meson sector 



13] , 14[ . As in the heavy meson sector on remains with several unsettled questions: 



1. Are there any special reasons for such big QCD "corrections"? 

2. Can we trust the QCD sum rule predictions and the (^-expansion 
when the a^-corrections are so big? 

3. How big are the a|-corrections? Is it possible to estimate them? 
These questions should be clarified in the near future. 



4.1 Residues and QCD sum rules 

To proceed with the usual QCD sum rules analysis, we evaluate the scalar correlator 
function P(uj) using the theoretical result Pope(^) given in Eq. (Q) and equate this to the 
dispersion integral over the contributions of hadron states. These consist of the lowest 
lying ground state with bound state energy A plus the excited states and the continuum. 
To leading order in 1/mq the bound state energy of the ground state is defined by 

^baryon = TUq + A, (24) 

where mq is the pole mass of the heavy quark. Note that the leading order sum rules 
do not depend on uiq at all since the heavy mass dependence has been eliminated by 
employing the heavy mass expansion. 

We assume that the continuum is given by the OPE expression above a certain thresh- 
old energy E c 0. For the hadron-side (h.s.) contribution to the spectral density we thus 
write 

Ph.s. M = Pg.s. (CO) + Pcont M , (25) 

where the contribution of the lowest-lying ground state (g.s.) baryon is contained in p g s . 
and is given by 

Pg . s {u) = l -FH{u-K). (26) 

In this expression F is the absolute value of one of the residues Fi (i = A, S, £*) of the 
baryonic currents defined by 

(0|J|A Q ) = F a u, (0|J|S o ) = F sM and (0| J V \E* Q ) = A^F^u v , (27) 

where u and are the usual spin- 1/2 and spin-3/2 spinors. Note that F-^* coincides with 
in the lowest order of the heavy quark mass expansion that we are working in. 
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As is usual we assume hadron-parton duality for the contribution of excited states and 
continuum contributions and take p CO nt(^) = 6{w — Ec)p(oo), where p is the result of the 
OPE calculations given in Eqs. (f|) and (|16|) . With these assumptions we arrive at the 
sum rule 



)F 2 r°° p(uj')duj' 

A — uj — iO Je c uj' — lu — iO 



PopeM= a ^ - +/ r ' n ( 2 8) 



or 

A-cj-iO Jo u'-u-iO p v ; ' v ; 

where the power counting part P pc (uj) is defined as the Fourier transform of that part 
of the correlator function P(t) which contains non-negative powers (t 2 ) n (n > 0). Finally 
we apply the Borel transformation 

uj n ( d\ n 

Bt = lim — — n, — uo — > oo (T = —uo/n fixed) (30) 

r(n) \ duj I 



to the sum rule in Eq. fl29|) . Using Bt(1/(u> — oj')) = exp(— lu'/T)/T we obtain the Borel 
sum rule 

l F 2 {p)e- K ' T = f EC p(u;',p)e-"'/ T du' + BP p . c XT) =: K(E c ,T,p), (31) 
I Jo 

where we reintroduced the //-dependence of the spectral density, which causes a p- 
dependence for the residue. The Borel-transformed BP P _ C ,(T) can be obtained directly 
from Pp. c .(t) by the substitution t — * — i/T (see the discussion in |L7]])- Note that the 
bound state energy A can be obtained from the sum rule in Eq. (|3~TD by taking the loga- 
rithmic derivative with respect to the inverse Borel parameter according to 

A = d\n(K(E c ,T,p)) 
Returning to the sum rule in Eq. ([IT]), one has 



l -F\p)e-^ = ™ 



(33) 



T6 ( h{xc) + S (( ln - f ^ Xc) ) ri + r2 )) 

+cE i G T 2 f 1 (x c ) + E 6 Q exp (-- ^ 

with the polynomials ri and r 2 presented in Eq. (|20| ) and the functions 

/ n (x) := / —e- x 'dx' = 1 - e~ x ^ — , 
Jo n\ ^ m\ 

f ( x ) ■= j X °^-]nx'e- x 'dx'. (34) 
Jo n\ 

In order to simplify the notation we have introduced the abbreviations 

Ic:= ^, £o:= - (^-^ and == ^f^. (35) 

The numerical analysis of the Borel sum rule is the subject of the section 5. 
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4.2 Anomalous dimensions 



The one-loop renormalization of the effective heavy baryon currents was considered in |T7] 



the two-loop case was studied in p4| . In general they differ from those in conventional 
QCD. The one-loop anomalous dimension of baryonic currents, namely the first coefficient 
in the expansion 7 = Yjk( a s I ^) k lki only depends on n and is given by JIT], 



7l = -l((n-2) 2 + 2). (36) 
The general (n, s)-dependent formula for the two-loop anomalous dimension case is 



rather lengthy and can be found in p4[ . As an illustration we list explicit values for the 
two-loop anomalous dimensions as calculated in the MS-scheme using the naive 75-scheme. 
One has (with explicit values given for Nf = 3) 



7ai 



7A2 



7si 



" 8 (S) + ^ (16C(2)+40iV/ - 796) O 2 ' (37) 



i-72.19 



- 4 (S) + ^ (16C(2)+20iV/ - 322) (S) 2 ' (38) 

ft:-26.19 

~ 4 © + l m{2) + 20Nf - 290) (S) 2 ' (39) 



i-22.63 



2 



7S2 = -^('^+^(48C(2) + 8JV> + 324)(^Y . (10) 



3 \4ttJ 27 v ' w J ' \4nJ 

" v ' 

«15.81 

4.3 Renormalization group invariant sum rules 

It is clear that the currents </(/i) depend on the renormalization scale /i. This dependence 
can be expressed by the renormalization group equation 

, d . w , dlnZj , , 

0*3- +7^=0. ^ := ^T— "> ( 41 ) 
a\i dm/i 

arising from the scale independence of the bare current J = Zj(n)J(n), where 7 is the 
anomalous dimension of the current discussed in the preceding subsection. To construct a 
renormalization group invariant quantity J- mv , the renormalized current J(n) is multiplied 
by some Wilson coefficient C(as(/i)), J\ m = J(fi)C(as(n)), which is subject to the "dual" 
renormalization group equation (see also (TTJ], [TBI ) 

(ji^--'y)C{as(ji)) = Q MM^ - 7K))C(a s ) = 0, (42) 

where /3 := d In a^/d In /x = Y,k( a s/^) k ^k is the beta function of QCD with 

2 19 
ft = -2(11 - -N f ) and ft = -4(51 - —N f ). (43) 
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The formal solution of Eq. (021) is given by 



C(a sM )=exp^ -^j. (44) 

Finally, the perturbative expansion of the beta function and the anomalous dimension up 
to second order in % gives 



4tt A V 71 (3 



(45) 



The first factor in Eq. ( [35| ) is the result of resumming the leading logarithmic terms 
(agln/i)™, where the result is valid only in the logarithmic approximation. As Eq. (|45|) 
shows, one needs to know also the two-loop anomalous dimension of the baryon current in 
order to obtain the evolution at next-to-leading log accuracy, e.g. in the order as{as hi ^) n - 
The usage of the invariance property of J; nv also provides a connection between currents 
at different renormalization scales, 

J(li 2 )C(a s (p 2 )) = J{ni)C{a s {iii)) => 
J(/i 2 ) = J{^i)C(a s (^i))C(a s (^ 2 )y 1 =: J(ni)U(fj,i,(j, 2 ), (46) 
r a s(Mi) da 7(a) \ 



U(ni,n 2 ) = exp 



as0 2 ) a /5(a)/ 
' ggOfO ^ fi - "5(^2) 7i (l2 U7) 



where U(fii,fi 2 ) is perturbatively evaluated up to next-to-leading order in as (see also 
the discussion in 0, ^5], [26|j). 

As being evident from Eq. (p7|), also the residues are functions of the renormalization 
scale parameter /z, the functional form for this dependence is the same as for the currents. 
So one can construct the renormalization group invariant F inv = F(/j,)C(as(n)) by means 
of the same Wilson coefficient. A renormalization group invariant sum rule can then be 
constructed by considering the expression 

~if nv exp(-A/T) = K(E c ,T,n)C(a s (fi)) 2 =: K inv (E c ,T). (48) 

The theoretical part of the sum rule K(Eq,T, /1) depends on the renormalization scale 
\i through the QCD perturbative corrections which contain the logarithmic factor ln(/i). 
On the other hand, the left hand side of Eq. is independent of the renormalization 
scale p by construction, and thus the right hand side must also be renormalization scale 
independent. It is easy to check this to first order in as by introducing a second scale \J 
and writing 

«^ = l_£drt AM 4). (49) 
as\n) o7r 

Remembering that p(u, /z) in Eq. ( p0|) appears as an integrand of K(E C , T, /x), one obtains 
cancellations (in first order of as) of the logarithmic factors ln(/x) in p(ui, fi)C(as(p)) 2 and 
thereby in K(Ec,T, fi)C{as{p)) 2 ■ The cancellation occurs because of r\ = —271. In this 
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paper we will make no usage of the renormalization group invariant sum rule in Eq. (f£8p. 
Instead of this we analyse the sum rule (|3~ID at some fixed point /i' = 1 GeV in order to 
estimate the bound state energy A and the residue F(fi'). The value of the residue F(fi) 
at other scales can then be obtained by using the evolution function U(fi',fi), while the 
//-independent function F inv can be immediately obtained by multiplying with C(as(ti))- 

5 Numerical results 

Let us discuss the sum rule analysis in some detail. We start by discussing the sum 
rules without radiative corrections and execute the analysis in consecutive steps. As 
Eq. (|33|) shows, the analysis of the non-radiatively corrected sum rules does not depend 
on which of the two different current cases are being discussed. First, we analyse the 
dependence of the bound state energy A on the threshold parameter Eq and the Borel 
parameter T in a large window of parameter space. The aim is to try and find regions 
of stability in T and Ec- By looking at the three-dimensional plots for A as functions 
of T and Ec we found a stability of the sum rules only in the case of the exponential 
ansatz for the non-local operator (q(0)q(x)). Keeping in mind the rather narrow window 
of 0.3 GeV < T < 0.4 GeV mentioned in the connection with Eq. (||) for the consecutive 
replacement r — > it — ► 1 /T, one ends up with a more reasonable discussion of the stability. 
We mention that the range of acceptable values for T is extended down to T > 0.2 GeV 
when radiative corrections are included, which enlarge the perturbative contributions. 
This, however, does not bring in a new region of stability. 

Returning to the analysis of the non-radiatively corrected sum rules for the Ag-baryon, 
we find areas of stability around Ec = 1.2 GeV in the window 0.3 GeV < T < 0.4 GeV. 
The range of confidence for Ec is 1.0 GeV < Ec < 1.4 GeV. Therefore in Fig. 2(a) 
we show plots for five values of Ec around E c est = 1.2 GeV, namely for Ec = E c est , 
Ec = E c est ±0.1 GeV and E c est ±0.2 GeV. From these curves we then can read off values 
for Ec and A with good sum rule stability, namely 

A(A) = 0.78 ± 0.05 GeV in the range E C (A) = 1.2 ± 0.1 GeV, (50) 

where the quoted errors present rough error estimates taken from Fig. 2(a) according to 
the interval for Ec in Eq. (|50|). 

Next we estimate the value of the residue. The sum rules depends on the three 
parameters A, E c and T. In Fig. 2(b) we plot \F\\ for a fixed bound state energy 
A(A) = 0.78 GeV in the indicated window for the Borel parameter T. The five different 
curves again correspond to the above five different values of Ec- Sum rule stability is 
found at 

|F A | = 0.023 ± 0.002 GeV 3 , (51) 

where the errors again represent rough error estimates taken from Fig. 2(b). 

Next we take into account the as-correction to the spectral density. As is evident 
from Eq. (pH|), the sum rule analysis now depends on which of the two types of baryonic 
currents are used. The results for the bound state energy for both cases are displayed in 
Fig. 2(c). Using the same analysis as for Fig. 2(a) we obtain 

A (A) = 0.78 ± 0.05 GeV in the range E C (A) = 1.1 ± 0.1 GeV. (52) 
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[17| 





@ 


L.O. 


N.L.O. 


E C (A) 


1.20 


1.20 ±0.15 


1.2 ±0.1 


1.2 ±0.1 


1.1 ±0.1 




1.46 


1.30 ±0.15 


1.4 ±0.1 


1.3 ±0.1 


1.3 ±0.1 


A(A) 


0.78 


0.79 ±0.05 


0.9 ±0.1 


0.78 ±0.05 


0.78 ±0.05 


A(E) 


0.99 


0.96 ±0.05 




0.90 ±0.05 


0.95 ±0.05 


A(£)-A(A) 


0.21 


0.17 




0.12 


0.17 


I*a| 


2.3 ±0.5 


1.7 ±0.6 


2.5 ±0.5 


2.3 ±0.1 


2.8 ±0.2 


1^1 


3.5 ±0.6 


4.1 ±0.6 


4.0 ±0.5 


2.6 ±0.2 


3.9 ±0.3 



Table 2: Sum rule results on non-perturbative and sum rule parameters of heavy ground 
state baryons. The continuum threshold parameter Eq, the bound state energy A and the 
difference between the two bound state energies are given in GeV, whereas the residues 
are listed in units of 10 -2 GeV 3 . The value of the Borel parameter is T = 0.35 GeV. 



Here we note the nice technical effect that the as-corrected sum rule is more stable at 
the lower value E c = 1.1 GeV for the continuum threshold but predicts the same bound 
state energy A. So there occurs some "stabilization" in A. 

Using the central value A(A) = 0.78 GeV one can then obtain values for the residue 
looking at Fig. 2(d), which give rise to the value 

| F A | = 0.028 ± 0.002 GeV 3 . (53) 

Doing the same analysis for the E baryon, at leading order in as we obtain 

A(E) = 0.90 ± 0.05 GeV, F C (S) = 1.3 ± 0.1 GeV and 
|F S | = 0.026 ± 0.002 GeV 3 (54) 

and including the as radiative corrections we have 

A(E) = 0.95 ± 0.05 GeV, F c (£) = 1.3 ± 0.1 GeV and 
|F S | = 0.039 ± 0.003 GeV 3 (55) 

The results are displayed graphically in Fig. 3(a,b) and in Fig. 3(c,d), respectively. 

Our predictions for the bound state energy A combined with the experimental charm 
and bottom baryon masses may be taken to calculate the charm and bottom quark pole 
masses rriQ. Taking into account the experimental results as given by the Particle Data 
Group |27|, namely ra(A c ) = 2284.9 ± 0.6 MeV, m(S+) = 2453.5 ± 0.9 MeV and m(A&) = 
5641 ± 50 MeV, we obtain the pole masses m c « 1500 MeV and rrib ~ 4860 MeV for the 
heavy quarks. The experimental difference of m(A c ) — m(S c ) w 167 MeV E7[ is quite 



near to our prediction m(A) — m(S) m 170 MeV. Here we present only central values. 
As was discussed above, the accuracy of our predictions is connected with the internal 
accuracy of the QCD sum rules method (mainly because of the dependence on the energy 
of continuum) and is probably not better than 20%. 

All the results are summarized in Table 2 where we compare our results with the 



leading order results obtained in [17, 19, 28|. This concludes our analysis. 
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6 Conclusions 



We have considered the Operator Product Expansion of the correlator of two static heavy 
baryon currents at small Euclidian distances and determined the a,s radiative corrections 
to the first Wilson coefficient in the expansion. Based on this expansion we formulated 
and analyzed heavy baryon sum rules for the A-type and E-type heavy baryons using two 
different types of interpolating fields for the baryons in each case. We have discussed in 
some detail the scale independence of the as sum rules which requires the consideration 
of the anomalous dimensions of the heavy baryon currents at the two-loop level. 

Similar to the case of heavy mesons the QCD radiative correction to the first term in 
the OPE is quite large and amounts to a 100% change in the perturbative contribution. 
The radiative correction to the perturbative term increase the calculated sum rule values 
for the baryon masses by about 10% and the residues by about 20 — 50% relative to the 
corresponding lowest order values. The sum rule results do not depend very much on 
which of the two possible interpolating fields is used in each case. The sum rule analysis 
is, however, quite sensitive to changes in the assumed threshold energy of the continuum. 
This sensitivity is the main source of uncertainty in our results and is partly due to the use 
of diagonal correlators. QCD sum rules based on the diagonal correlators feature a leading 
order spectral density which grows rapidly as p{u) ~ uj 5 . This rapid growth introduces a 
strong dependence of the sum rule results on the assumed energy of continuum. Second, 
the QCD radiative correction to the leading order spectral density is about 100% at the 
renormalization scale /i — 1 GeV. We may try to make the coefficient at 0:5 in these 
corrections to be moderate considering the very low renormalization scale /i — 10 MeV . 

We have not considered non-diagonal sum rules which come in when one considers 
correlators between two different currents with the same quantum numbers. These non- 
diagonal sum rules bring in some new features such as a more "normal" behaviour of the 
spectral density p{uj) ~ (qq)oo 2 and thus probably more moderate QCD corrections to 
this spectral density. On the other hand, the leading term for non-diagonal sum rules is 
proportional to the quark condensate, whose value (qq) = (—0.23 ± 0.02 GeV) 3 is known 
only with an accuracy of 10%, which gives an additional uncertainty in the result for 
non-diagonal sum rules. The analysis of the non-diagonal sum rules will form the subject 
of a subsequent paper. 
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Appendix A: Diagrammatic contributions 



In this Appendix we collect together results on the calculation of the two-loop and three- 
loop contributions to the correlator of two heavy baryon currents. We start with the 
two-loop contribution depicted in Fig. 1 (i — 0) where one has 

*>« = (C-4KJ-3) - (A1) 

We have introduced the abbreviation E n = r(l — e)"T(l + ne) (with natural numbers 
n — 1, 2, 3, ... ) which is also used in the subsequent presentation of the three-loop results. 
For the three- loop contributions i — 1,2 and 4 depicted in Fig. 1 one has 

2(D -2){2D - 7)E 3 ~ - 

1 ~ 9(3D - 11) (3D - 10) (3D - 8) (3D - 7)E 2 1 Tr ^ Tt > 



with 



~ 2{D - 2)E\ (D-2)(3P-10)£ 2 

1 ' (D-4) 3 (D-3) 2 (D-4) 3 (D-3) 2 (2D-7)' 1 ' 



(D - 2)E 2 



2(D-4) 2 (D-3)(2D-7)' 

-(D-2)E 2 
(D-4) 2 (D-3) 2 (2D-7)' 



The contribution of diagram (2) in Fig. 1 is the most involved one. In order to be able to 
write the results in a compact form we introduce the abbreviations 

Q 1 = r(l -e) 2 r(l + e)/r(l -2e) and 

Q 2 = r(l-e) 3 r(l + 2e)/r(l-3e). (A3) 
In terms of the basic structure terms 

f = Tr(IW), f x = T^f 7^7") and f 2 = Tr(f 7^7,^7 "7"), (A4) 
one obtains 

B2 = 36(D-3)(3D- 3 ll)(3D-7)g 2 g^ With 



^2,0 



12(D-2) 2 g 2 24D(D - 2) 2 Q 2 



(D - 4) 3 (D - 3) 2 (D - 1) (D — 4) 3 (D - 1)(3D - 10)(3D - 8) ' 



~ (D 2 -7D + 16)Q 2 4 (D 2 -4D + 8)g 2 

2,1 (D-4) 2 (D-3) 2 (D- 1) (D-4) 2 (D- 1)(3D- 10)(3D-8)' 1 ' 



2,2 



SQ( 4Q 



2 



(D - 4) 2 (D - 3)(D - 1) (D - 4)(D - 1)(3D - 10)(3D - 8) ' 
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Figure Captions 

Fig. 1: Two- loop and three-loop contributions to the correlator of two heavy baryon 
currents. (0) two-loop lowest order contribution, (l)-(4) three-loop 0(0:5) 
contributions. 

Fig. 2: Sum rule results on non-perturbative parameters of the Aq as functions of 
the Borel parameter T. Shown are five equidistant curves centered around 
E c = E c cst with a distance of 100 MeV, E c growing when going from 
bottom to top. These are in detail 

(a) lowest order sum rule results for the bound state energy A(A) 

(b) lowest order sum rule results for the absolute value of the residue F\ 

(c) 0(as) sum rule results for the bound state energy A(A) 
for currents Jai (solid) and J\ 2 (dashed) 

(d) 0(a s ) sum rule results for the absolute value of the residue F A 
for currents Jai (solid) and J A2 (dashed) 

Fig. 3: Sum rule results on non-perturbative parameters of the Sq as functions of 
the Borel parameter T. Shown are five equidistant curves centered around 
Ec = E c cst with a distance of 100 MeV, Ec growing when going from 
bottom to top. These are in detail 

(a) lowest order sum rule results for the bound state energy A(S) 

(b) lowest order sum rule results for the absolute value of the residue F% 

(c) 0(a,s) sum rule results for the bound state energy A(E) 

for the current doublets {Jsi, Js*i} (solid) and { Js2, ^2*2} (dashed) 

(d) 0(a s ) sum rule results for the absolute value of the residue F s 
for the corrent doublets {Jsi, Je*i} (solid) and { J S2 , Js*2} (dashed) 
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